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Abstract 

Investigation of long-term meteorological satellite data revealed statistically significant 
vegetation response to climate drivers of temperature, precipitation and solar radiation with 
exclusion of fire disturbance in Alaska. Abiotic trends were correlated to satellite remote 
sensing observations of normalized difference vegetation index to understand biophysical 
processes that could impact ecosystem carbon storage. Warming resulted in disparate 
trajectories for vegetation growth due to precipitation and photosynthetically active radiation 
variation. Interior spruce forest low lands in late summer through winter had precipitation deficit 
which resulted in extensive fire disturbance and browning of undisturbed vegetation with 
reduced post-fire recovery while Northern slope moist alpine tundra had increased production 
due to warmer- wetter conditions during the late 1990s and early 2000s. Coupled investigation of 
Alaska’s vegetation response to warming climate found spatially dynamic abiotic processes with 
vegetation browning not a result from increased fire disturbance. 
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23 1.0 Introduction 


24 The global warming imprint has left an indelible mark on Arctic terrestrial processes 

25 observable from satellite remote sensing instruments [Kimball et al., 2007]. Future decades hold 

26 high latitude changes from abiotic impacts to biophysical processes. Mean annual surface 

27 temperature increased > 2° C in Alaska during the past 25 years which has been observed to have 

28 marked regional impacts altering ecosystem functioning through processes of snow melt 

29 modifying albedo and/or hydrology [ Chapin et al., 2005; Dye and Tucker, 2003; Sturm et al., 

30 2001], thawed permafrost increasing depth of active layer [Osterkamp and Romanovsky, 1999], 

3 1 enhanced vegetation growth via growing season extension [Epstein et al., 2004; Jia et al., 2003; 

32 Stow et al., 2003; Walker et al., 2003], increased fire disturbance periodicity due to summer 

33 drought [Stocks et al., 2003], reduced growth from temperature induced drought stress [Barber et 

34 al., 2000], and changed shrubland cover [Tape et al., 2006]. Increased regional warming could 

35 produce carbon sinks [Sitch et al., 2007] or sources [Goetz et al., 2005; Pisaric et al., 2007] due 

36 to disturbance interval [Amiro et al., 2001] or variance in cycling rate [Kimball et al., 2007] 

37 which could feedback to global climate. 

38 Ecosystem dynamics altered by climate change are complex in cold regions, which has 

39 produced future carbon status uncertainty [Comelissen et al., 2007]. Warming in Alaska has 

40 accelerated from 0. 1 5 ± 0.02 to 0.3 ± 0. 1 per decade and now Arctic summers are warmer than 

41 400 years prior, resulting in large impacts to water dependent processes [Barber et al., 2004; 

42 Chapin et al., 2005; Riordan et al., 2005]. Alaska is an ideal to understand northern biomes as it 

43 exhibits increased Arctic slope vegetation productivity while negative trends exist in the interior 

44 [Neigh et al., 2008; Verbyla, 2008]. Climate and biophysical ecosystem interaction change is 

45 critical to understand as it may be indicative of future high-latitude processes. To address this 



46 the hypothesis herein was long-term summer climate warming has produced photosynthetic 

47 trends either positive or negative depending on regional temperature and precipitation 

48 requirements while fire and/or insect outbreak disturbance processes are second order. Prior 

49 studies have not provided spatially contiguous correlation results of abiotic forces on vegetation 

50 growth with fire disturbance regimes in Alaska. This study seeks to understand ecosystem 

51 dynamics observable through satellite measurements. 

52 2.0 Experiment Design and Data 

53 Multiple geospatial datasets were acquired to understand vegetation production variance 

54 throughout Alaska. Distinguishing between multiple processes and feedbacks is difficult 

55 considering they are often driven by one another, for example warming drought-stress inducing 

56 fire followed by insect-outbreak. To illustrate Alaskan vegetation processes, a figure was 

57 developed for disturbance agents to vegetation production indicated with boxes, and feedbacks 

58 shown as ovals (Figure 1). To further explore intra-seasonal dynamics, data trends were 

59 calculated on monthly values. 

60 Two experiments were conducted to understand climate-vegetation relationship and 

61 implication of fire disturbance. Experiment one calculated per-pixel correlation between 

62 normalized difference vegetation index (NDVI) the measure of photosynthetic capacity of 

63 vegetation, surface air temperature, precipitation, and photosynthetically active radiation data 

64 (PAR). Experiment two examined ecoregion mean of positive or negative significant abiotic - 

65 NDVI correlation with and without fire disturbance. 



66 2.1 Geospatial & Disturbance Data 

67 Historical burned area perimeters were derived from the Alaska Fire Service GIS Group, US 

68 Department of the Interior, Bureau of Land Management \AFS, 2008] 

69 [http://agdc.usgs.gov/data/blm/fire/index.html]. Burned area data quality varied to 1945 due to 

70 development method and resources available. A marked increase in burned area from the 1980s 

71 ~1 .6 million hectares to 2000s ~6.5 occurred throughout interior Alaska. This change in fire 

72 history could impact ecosystem productivity and/or recovery observed with coarse resolution 

73 AVHRR. Burned area was reprojected to Global Inventory Modeling and Mapping Studies 

74 (GIMMS) North America Albers projection and converted to fractional 64 km from 2 km using 

75 average pixel aggregations. 

76 Historical insect outbreak data were derived from Alaska GIS Group, US Department of the 

77 Interior, Bureau of Land Management [ USFS , 2008] [http://agdc.usgs.gov/data/projects/fhm/#K] and 

78 converted to 8 km average pixel aggregations. Outbreak area was collected for > 40 disturbance 

79 types with varying implications to vegetation health from reduction in leaf area inducing early 

80 senescence to mortality. Aspen leaf miner was the most extensive insect disturbance throughout 

81 interior and southeast regions during the 2000s although it does not cause mortality [ Wagner et 

82 al., 2007]. However it could impact NDVI measurements [ Verbyla , 2008]. All insect 

83 disturbances were prevalent in bottom lands of closed spruce hardwood forests, and open, low 

84 growing spruce forest neighboring river banks. 

85 To spatially understand remote observations of disturbance, ecoregions were subset based 

86 upon US Department of Agriculture Forest Service ECOMAP Version 2.0 [Nowacki et al., 

87 2001], to 34 sub-regions based upon dominant vegetation cover, climate, and altitude. 

88 Investigation sought to understand if disturbance processes initiated interior NDVI decline and 



89 ECOMAP provided a means to subset National Oceanic Atmospheric Administration’s polar 

90 orbiting (NOAA) satellite measurements of vegetation photosynthetic capacity. Most ecoregions 

91 include large samples > 130, 8-km pixels (Table 2). 

92 2.2 Remote Sensing and Climate Data 

93 GIMMS version ‘g’, 1 982 to 2005 bimonthly AVHRR NDVI data [ Tucker et al., 2005] 

94 provide a consistent inter-calibrated record for long-term vegetation studies. These data were 

95 corrected to account for orbital drift, minimize cloud cover, compensate for sensor degradation, 

96 and stratospheric volcanic aerosols effects [Tucker et al.]. GIMMS Alaska data contains nearly 

97 ~25,000 8 km pixels extending back to 1982 from 2005 yielding 144 growing season months for 

98 correlation. July NDVI > 0.5 was used as a threshold to exclude glaciers and sparsely vegetated 

99 mountainous regions from calculations. 

100 Monthly climate data were derived from Leemans & Cramer climatology [ Leemans and 

101 Cramer, 1991] and GISSTEMP anomalies [Hansen et al., 1999], solar radiation from the 

102 International Satellite Cloud Climatology Project (ISCCP) [Bishop and Rossow, 1991], and 

103 precipitation from the Global Precipitation Climatology Project version 2 (GPCP) [Adler et al., 

104 2003]. Climate data were detrended and reprojected to GIMMS North America Albers 

105 projection and bilinear interpolated to AVHRR NDVI grid cell resolution. Climate trends are 

106 calculated in a similar manner as prior NDVI investigations [Slayback et al., 2003] with a least 

1 07 squares linear fit per pixel from 1 982 through 2005 applied with pixels having a significance of 

1 08 less than 0.1 or a confidence of 90% retained. Values presented are slope multiplied by 24- 

1 09 years, between 1 982 - 2005. Correlation of climate to NDVI has been performed in numerous 

110 studies to understand implications of abiotic changes to ecosystems [Braswell et al., 1997; 



111 Myneni et ah, 1996; Neigh et al., 2007; Potter and Brooks, 1998]; similar methods are employed 

112 herein to understand browning of interior Alaska. 

113 3.0 Results & Discussion 

114 Calculation of monthly mean climate and trends correlated with NDVI from 1 982 through 

115 2005 revealed changes in productive growing season. North Slope mean July temperature 

116 ranged ~5 - 1 0 °C with prior and later months experiencing temperatures close to freezing 

117 potentially allowing snow cover altering growing season NDVI measurements. It is considered a 

118 polar desert with most precipitation occurring late growing season and seasonal snow thaw 

119 contributes to early productivity. Warmest Alaskan summer temperatures > 1 5 °C occur in the 

120 continental climate of the interior and monthly precipitation was low < ~20 mm throughout early 

121 and mid-growing season increasing to the southeast > ~150 mm in September. Available High- 

122 latitude PAR varies markedly in Alaska >150 W/m 2 in June, and < 50 W/m 2 in September. 

123 Note mean seasonal vegetation growth is greatest in July during temperature maxima, while 

124 having limited water availability throughout interior and northern Alaska. Abiotic driver change 

1 25 could alter early and mid-season production depending upon temperature impact on vapor 

126 pressure deficit. 

127 Alaska exhibited strong Arctic Slope vegetation growth with declining interior trends (Figure 

128 2). Temperature trends > 2 °C are prevalent during early months with little late season variation. 

1 29 All of Alaska experienced May - June increased precipitation totaling > ~20 mm from 1982- 

1 30 2005 except for the southeast handle which has the greatest precipitation. North Slope and 

131 southeastern coastal Alaska had increased late season precipitation, while a declining interior 

132 trend extended through winter months. PAR changes were minimal <±10 W/m . 



133 Experiment one found abiotic drivers of temperature, and precipitation had intra-seasonally 

1 34 moderate positive or negative relationship to vegetation productivity, while PAR had weak 

135 correlation. Temperature correlation to NDVI revealed > 0.5 in May and June throughout 

136 Alaska indicating earlier growing season start (Figure 3). Precipitation was positively correlated 

137 > 0.4 to vegetation growth in May and June on the North Slope, while late season precipitation 

138 decline and winter snowpack correlated > -0.5 negatively. 

139 Experiment two calculated NDVI correlation to abiotic variables in locations of fire 

140 disturbance > 50%, 8-km pixel burned (Figure 3, inset black bars) and with fire exclusion 0%, 8- 

141 km pixel burned (Figure 3, inset white bars). Subtle ecoregion difference between bumed- 

142 unbumed occurred early season with larger difference during mid to late season months. 

143 Negative NDVI interior regions had moderate negative correlation with surface temperature and 

144 precipitation, with weak positive correlation to PAR within unbumed sites and slightly stronger 

145 negative correlation in burned sites. Most ecoregions had similar correlation whether burned or 

146 unbumed with few ecoregions exhibiting stronger negative precipitation correlation and less 

147 positive correlation to temperature in burned sites. 

148 Regional vegetation growth and browning due to climate change occurred from 1982 through 

149 2005. Results presented are similar to prior North Slope reports with temperature increases 

1 50 driving growth in vegetation; however increased precipitation was also found to have an impact. 

151 Suspected interior lowland drying appeared during late season and through reduced snowpack. 

1 52 Satellite observations found spatially contiguous regions of vegetation productivity change 

1 53 which had moderate relationship to temperature and/or precipitation. North Slope vegetation 

1 54 growth appeared from warmer- wetter conditions while interior drying-browning vegetation 

155 appeared late growing season, followed by reduced winter snow pack leading to reduced spring 



156 moisture. Coupled warming-drying climate with poor post disturbance recovery is suspected 

1 57 driver for early season browning. Warming permafrost could contribute through nutrient cycling 

158 and surface saturation [ Chapin et al., 2005]. However, capturing active layer dynamics is 

1 59 beyond the scope of this investigation. No marked mid-growing season interior climate trends 

1 60 were found. Fire and insect outbreak could reduce mid-season vegetation productivity although 

161 correlation difference between disturbed and undisturbed sites was minimal. Increased fire 

1 62 disturbance interval appeared not to cause browning, but is a result of long-term drying. 

1 63 Correlated data revealed regional climate change could impact vegetation production and 

1 64 Alaskan terrestrial carbon cycle balance. Long-term spatial climate records appear to be robust 

1 65 using simple correlation significance to understand climate influence on vegetation growth. 

166 Future investigation will quantify regional carbon budget disturbances in ecosystem simulations. 
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Table 1. ECOMAP ecoregions (34) used in spatial analysis of 
abiotic variables and NDVI. Code numbers referred to in Fig 3 . 


Ecoregion 

Code 

# Pixels 

Area (10 3 km 2 ) 

Beaufort Coastal Plain 

l 

818 

52.4 

Brooks Range Foothills 

2 

1659 

106.2 

Brooks Range Mountains 

3 

2345 

150.0 

Olgivie Mountains 

4 

641 

41.0 

Ray Mountains 

5 

992 

63.5 

Yukon Flats 

6 

528 

33.8 

Upper Kobuk-Koyukuk 

7 

872 

55.8 

Kotzebue Sound Lowlands 

8 

338 

21.6 

Nulato Hills 

9 

815 

52.2 

Seward Mountains 

10 

778 

49.8 

Yukon Bottomlands 

11 

1013 

64.8 

Dawson Range 

12 

1021 

65.3 

Yukon-Kuskokwim Delta 

13 

1367 

87.5 

Kuskokwim Colluvial Plain 

14 

791 

50.6 

Kuskokwim Mountains 

15 

1142 

73.1 

Alaska Range 

16 

1431 

91.6 

Copper River Basin 

17 

270 

17.3 

Cook Inlet Lowlands 

18 

427 

27.3 

Wrangell Mountains 

19 

394 

25.2 

Nushagak-Lime Hills 

20 

558 

35.7 

Northern Chugach Range 

21 

227 

14.5 

Northern Aleutian Range 

22 

579 

37.1 

Chugach Range 

23 

295 

18.9 

Kenai Mountains 

24 

161 

10.3 

St. Elias Range 

25 

489 

31.3 

Ahklun Mountains 

26 

531 

34.0 

Gulf of Alaska Fjordlands 

27 

350 

22.4 

Boundary Range 

28 

311 

19.9 

Gulf of Alaska Forelands 

29 

133 

8.5 

Alexander Archipelago 

30 

870 

55.7 

Bristol Bay Lowlands 

31 

635 

40.6 

Alaska Peninsula 

32 

481 

30.8 

Aleutian Islands 

33 

193 

12.4 

Aleutian Colluvial Plain 

34 

179 

11.5 


Figure 1. Spatial-temporal vegetation productivity dynamics with interacting feedbacks in Alaska. 

Figure 2. (A) May, (B) June, (C) July, (D) August, and (E) September 1982-2005, slope multiplied 24-years of 
NDVI, surface air temperature, precipitation, and photosynthetically active radiation. 

Figure 3. ECOMAP ecoregions, AFS burned area perimeters converted to fraction of 8 km pixel displaying severe 
bum years 2004-2005 overlaid upon a digital elevation model. (Lower) Mean of positive or negative deseasonalized 
detrended significant correlation coefficients for the spatial regression of NDVI versus temperature (red bars), 
precipitation (blue bars), and photosynthetically active radiation (orange bars) by month and ecoregion presented as 
mean of entire ecoregion, without fire (white inset bars), and bum locations > 50% of 8 km 2 pixel (black inset bars). 
(A) May, (B) June, (C) July, (D) August, and (E) September 1982-2005, values with a significance > 0.05 and July 
NDVI < 0.5 excluded. 




Figure 1 . Hypothesized spatial-temporal vegetation productivity dynamics potentially observed with 
meteorological satellite measurements with interacting feedbacks in Alaska. 
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Figure 2. (A) May, (B) June, (C) July, (D) August, and (E) September 1982-2005, slope multiplied 24-years of 
NDVI, surface air temperature, precipitation, and photosynthetically active radiation. 
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Figure 3. ECOMAP ecoregions, AFS burned area perimeters converted to fraction of 8 km pixel displaying severe 
bum years 2004-2005 overlaid upon a digital elevation model. (Lower) Mean of positive or negative deseasonalized 
detrended significant correlation coefficients for the spatial regression of NDVI versus temperature (red bars), 
precipitation (blue bars), and photosynthetically active radiation (orange bars) by month and ecoregion presented as 
mean of entire ecoregion, without fire (white inset bars), and bum locations > 50% of 8 km 2 pixel (black inset bars). 
(A) May, (B) June, (C) July, (D) August, and (E) September 1982-2005, values with a significance > 0.05 and July 
NDVI < 0.5 excluded. 




